Low-lying continuum states of drip-line Oxygen isotopes 
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Low-lying continuum states of exotic oxygen isotopes are studied, by introducing the Continuum- 
Coupled Shell Model (CCSM) characterized by an infinite wall placed very far and by an interaction 
for continuum coupling constructed in a close relation to realistic shell-model interaction. Neutron 
emission spectra from exotic oxygen isotopes are calculated by the doorway-state approach in heavy- 
ion multi-nucleon transfer reactions. The results agree with experiment remarkably well, as an 
evidence that the continuum effects are stronger than ~1 MeV, consistently with the shell evolution 
in exotic nuclei. The results by this CCSM doorway-state approach are compared with calculations 
on neutron-scattering resonance peaks made within the CCSM phase-shift approach and also with 
those obtained in the Gamow shell model, by taking the same interaction. Remarkable similarities 
in peak energies and certain differences in widths are then obtained. 

PACS numbers: 21.60.Cs,24.10.Cn,24.30.Gd,27.30.-|-t 



I. INTRODUCTION 

Exotic nuclei far from the /3-stability line provide us 
with new interesting features. A good example is the 
neutron halo and other continuum-related properties of 
drip-line nuclei Another example can be the evo- 
lution of shell structure, including the change of magic 
numbers, due to characteristic properties of the nuclear 
force 0, as being increasingly observed 0, [Hi- Al- 
though there are many bound nuclei between /3-stability 
and drip lines showing the shell evolution, one eventu- 
ally approaches the drip line, by adding more neutrons 
(or protons). Low- lying states are then in the continuum 
even if the ground state is still bound. Thus, the physics 
of continuum and that of shell evolution should meet. 
It is of much importance to clarify how the coupling to 
the continuum affects various structures generated by the 
nuclear force. We shall study, in this paper, continuum 
spectra of neutron-rich O isotopes. The O isotopes show 
unusual shell structure and magic numbers (iV=14, 16) 
[1, IB-[§1, ending up with unbound 0^3/2 orbit. 

The neutron drip line is quite different between O and 
F isotopes which can be explained by unbound 

(bound) neutron 0^3/2 orbit for O (F) isotopes. However, 
despite their importance, properties of unbound neutron 
d3/2 states in O isotopes remain to be clarified. 

We extend the normal shell model (SM) so as to in- 
clude continuum states as a part of SM basis. By this 
continuum-coupled shell-model (CCSM) calculation, we 
can see precisely the location and shape of the peaks 
of various continuum spectra. Although the use of an 
infinite wall belongs to standard techniques, the wall 
is placed at unusually far distance so that the SM ba- 



sis includes many (^3000) discretized continuum states, 
both resonant and non-resonant. This allows us to pro- 
duce spectrum shapes accurately and treat the asymp- 
totic behavior of wave functions appropriately. A finite- 
range nucleon-nuclcon (NN) force, which is consistent 
with successful SM interaction for bound states, is intro- 
duced for the evaluation of continuum effects, resulting in 
significant energy shifts from the bound(-state) approxi- 
mation. Such novel approach leads us to quite interesting 
results, in agreement with recent experimental data. 



II. CONTINUUM-COUPLED SHELL MODEL 

As the 0^3/2 orbit remains unbound in oxygen isotopes, 
we generate discretized continuum single-particle basis 
states of ^3/2, by introducing an infinite wall on top of the 
Woods-Saxon potential. The single-particle Hamiltonian 
is then written as. 



Ho = T + Uws + V; 
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where T denotes kinetic energy and Uws is the Woods- 
Saxon potential with radius i?=1.09A^/^ fm, diffuseness 
0=0.67 fm. The depth parameter Vq will be discussed 
later. The standard LS term is included Here, T4,aU 
is an infinite wall placed at 3000 fm from the center of 
the nucleus so as to obtain a sufficiently high density of 
discretized states. As shown later, the final results do 
not depend on the position of the infinite wall if the wall 
is far enough. The i-th ^3/2 state is denoted as \id3/2) 

with i=l, 2, - • • ,iinax- 

Six neutrons are assumed, for simplicity, to be in the 
0^5/2 orbit. Namely, the inert ^^O core is assumed in the 
following CCSM calculations. Thus, the basis states of 
valence neutrons are constructed from the bound lsi/2 
state and the discretized continuum 1^3/2 (i = 1,2,...) 
states with various sets of occupation numbers. 
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The CCSM is proposed to treat correlations due to nu- 
clear forces in continuum states in a manner consistent 
with the SM for bound states. Realistic SM interactions 
are given usually in the form of two-body matrix ele- 
ments (TBME's). In order to calculate matrix elements 
of the Hamiltonian for the continuum basis states dis- 
cussed above, the SM interaction is extended to an an- 
alytic function, V, depending on variables of two inter- 
acting nucleons. As the realistic SM interaction contains 
many effects, from free-nucleon scattering to in-medium 
polarization corrections, it is virtually impossible to ex- 
press in a simple functional form for all single-particle 
orbits. On the other hand, the continuum coupling is 
important for few selected orbits in many cases. For 
instance, in the present case, neutron active orbits are 
lsi/2 and 0^3/2, and c?3/2 is extended to continuum basis 
states. We introduce a simple as a modeling of a realis- 
tic interaction within this small model space. Continuum 
*<^3/2 states differ among themselves in their radial wave 
functions, and are mixed due to the coupling by nuclear 
force, producing correlated eigenstates even in the con- 
tinuum. The coupling occurs primarily through the inner 
(i.e., resonance- like) part of such wave functions, because 
in the outer part nucleons are far from each other. The 
inner part resembles, more or less, bound 0^3/2; apart 
from the overall amplitude. The interaction with a neu- 
tron in lsi/2 is important too. Therefore, it is the mini- 
mal condition that V should reproduce relevant TBME's 
of a realistic SM interaction involving bound (i.e.. Har- 
monic Oscillator (HO)) lsi/2 and 0c?3/2 orbits. In this 
work, V is fixed by this minimal condition. 

As a simple yet general form oi V, a. superposition of 
two Gaussian functions with spin dependences is taken: 

E + (2) 

1=1,2 

where r is the inter-nucleon distance, and a implies spin 
of a nucleon. Here, di_2=1.4, 0.7 fm are prefixed for sim- 
plicity. The other parameters are fixed later. The CCSM 
Hamiltonian is then written as 

H = Ho + V = Y^ ijUj + V, (3) 
j 

where ij and Uj are, respectively, the single-particle 
energy (SPE) and the occupation number operator, 
and j denotes single-particle states including discretized 
continuum ones. 



III. CONSTRUCTION OF CCSM 
HAMILTONIAN 

We define the SM Hamiltonian to be used for fixing V. 
The SDPF-M Hamiltonian [l^l, a revision of USD is 
taken. The SDPF-M Hamiltonian has been shown to be 
successful for describing bound srf-shcU nuclei including 



those with a large number of neutrons. We begin with 
the fine tuning of the SDPF-M so that it can reproduce 
one neutron separation energy (S'„) of ^'^O exactly, by 
the following minor change of the monopole interaction, 

5(lsi/20d5/2|V^|lsi/20d5/2)^=' = -0.03 MeV. (4) 

Note that this change is less than 5% of the original value. 
While this fine-tuned Hamiltonian is for the full SM di- 
agonalization, we introduce a simpler SM Hamiltonian 
to be used in the filling configuration scheme. The filling 
configuration means that valence nucleons are put into 
the lowest possible orbit one by one, and the diagonal- 
ization is carried out within such configurations. This 
truncation makes sense if the valence orbits are well sep- 
arated. For oxygen isotopes with N > 14, this is the case 
as represented by N=14 and 16 gaps ^4 MeV [l7l |. 

The SPE's and TBME's appropriate for the filling 
scheme are determined so that the results of the fully 
mixed SM calculations are reproduced by the calculations 
in the filling scheme, as far as low-lying energy levels of 
23-26Q j-QioXive to the ^^O ground state are concerned. 
Note that the active model space is (lsi/2-0d3/2) for the 
filling scheme. The resultant Hamiltonian is referred to 
as modified SDPF-M. 

Next, the parameters of V are determined from the 
modified SDPF-M Hamiltonian. Relevant TBME's 

(lSl/20d3/2|l^|lSi/20(i3/2)j=l,2 (5) 
(Od3/20d3/2|-t^|0d3/20d3/2)j=0.2 (6) 

are calculated for each term of V by using HO single- 
particle wave functions. The parameters of V are op- 
timized by the fit to the corresponding TBME's of 
the modified SDPF-M, under the condition \V{r)\ < 
400 MeV. One can uniquely determine the four parame- 
ters, ai,gi (i = 1,2), by four constraints. If we do this, 
however, the parameters appear to be unnaturally large, 
while the final SM result turned out not to change much 
from the above one. We therefore choose the parame- 
ters fixed by the the fit mentioned above. The ob- 
tained values are 5i=39.9 McV, .g2=-409 MeV, ai=0.38 
and 02 =0.09. Their validity will be examined from a 
different viewpoint later. 

The Woods-Saxon depth parameter of Hq in cq. ([T]) is 
determined so that when the HO wave function is used 
for 0^3/2, the expectation value of Hq becomes equal to 
the 0^3/2 SPE of the modified SDPF-M Hamihonian. 
The same Woods-Saxon potential is used in all CCSM 
calculations in this paper as well as the other parameters 
in the Hamiltonian in eq. ([3]). The CCSM Hamiltonian 
is thus fixed. 



IV. CONTINUUM STATES IN ^^o 

We move to the structure of ^^O. As the filling con- 
figuration is taken, the ground state is composed of two 
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neutrons in the lsi/2 orbit. This orbit is well bound with 
S'„=2.7 MeV for ^^O [H], and hence is represented by a 
HO wave function. The excited state of the angular mo- 
mentum J (=1,2) lies in the continuum, and is given by 
a superposition of basis states. 



|iJ+) = |lSi/2 (8)1^3/2;^"*"), 



(7) 



where denotes the i-th discretized continuum state 
introduced earlier. In the present CCSM calculation, 
*max for corrcspouds to excitation energy of ^50 

MeV above threshold. After the diagonalization, we ob- 
tain the fc-th eigenstate (eigen energy Ek) of the angular 
momentum J, 



|^J+), 



(8) 



where c\ ' denotes mixing amplitude. 

We discuss a reaction leading to '^'^O continuum states 
through doorway state created by the sudden removal of 
a proton and a neutron from ^^F. We assume that after 
the removal of the proton, there is no proton in the sd 
shell in the doorway state. For ^^F (A^=17), six neutrons 
are assumed to remain in Qd^/2 of the ^^O core before 
and after the reaction, because of a lar ge g ap 4 MeV) 
between 0^5/2 and lsi/2 around N=1A [l7[. The remain- 
ing three neutrons are assumed to be in the lsi/2-0(i3/2 
space in the initial state. After the removal of one of 
the three neutrons, there are configurations: (a) 15^/2' 
(b) ls\^2^'^\/2^ ^'^'^ ('^) ^'^\/2- C*^)' ^^"^ initial state 
should be of either configurations Qd'l^^^ or Isy^Ody^i 
which are excluded in the filling scheme and should be 
indeed negligible because of a large lsi/2-0(i3/2 neutron 
gap in nuclei being studied 0, • We do not discuss the 
configuration (a), because this is nothing but the bound 
ground state of ^''O. Thus, only the configuration (b) is 
considered for the doorway state. Since ^^F is bound, 
lsi/2 and 0(^3/2 orbits of the doorway state are expressed 
by the HO wave functions. The lsi/2 wave function keeps 
this feature after the reaction, because it is well bound as 
confirmed experimentally by Kanungo et al. How- 
ever, the 0c?3/2 is not bound in ^^O, and the neutron 



^3/2 

0d}^^2 HO wave function becomes the initial state 



(doorway state) for the following time evolution. After 
the reaction, this two-ncutron system of ^^O is described 
by the CCSM Hamiltonian in cq. ([3]). Its eigenstate with 
in cq. dHl). 

with amplitudes 



the energy is | J^t) 

The isy2 ^^^3/2 doorway state is expanded by \Jjt) 



(J+|lSi/20d3/2;J-' 



and its overlap probability is given by. 



pi'' 



(^^3/210^3/2) , 

(9) 
(10) 
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FIG. 1: (color online) Probability of emission of ^3/2 neutron 
from 23-25Q ^ function of the energy of emitted neutron. 
Vertical bars, solid and dashed curves denote, respectively, the 
results of the bound approximation, the CCSM calculations 
with and without the residual interaction {V in eq. (|2])). 



Namely, with this probability, the doorway state evolves 
in time after the reaction as the state iJjt)- The \Jit) 
state represents a freely moving neutron at far distance 
with the kinetic energy Ek-Eo with Eq being the ^'^O 
ground-state energy. This is nothing but the emission of 
^3/2 neutron at the energy Ei^-Eq in the channel J+. Al- 
though the reaction phase space should be considered for 
the neutron emission cross section and it varies in prin- 
ciple as a function of Ek , this change is rather minor and 
is neglected, because of high bombarding energy. Thus, 
the probability p^f'' is considered to represent the spec- 
trum of d3/2 neutron emitted in the channel J"*" except 
for the overall absolute magnitude. We do not discuss 
the absolute magnitude of the cross section. 

Figure [1] shows ^^"''''8. While p^ff^ is defined for a dis- 
crete value of Ek originally, we smear it out by Gaussian 
function under the usual condition that the value of p^"^-* 
is equal to the integral of this Gaussian. The width of 
the Gaussian is varied in accordance with the interval of 
discrete energies E^s. Figures [T] (a,b) show the values of 
p'lf'' for J^=\ and 2, as functions of the energy of emit- 
ted neutron. The solid curve shows the present CCSM 
calculation. This curve should be compared to the ex- 
perimentally observed spectrum, except for the overall 
magnitude. 

Regarding the comparison to experiments, we first 
mention the absence of bound excited states in ^^O, con- 
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FIG. 2: (color online) Peak energies of one-neutron ('^"'"^^O) 
or two-neutron (^^O) emissions compared with experiment 
(Exp.). Calculations are made by bound states for SDPF-M, 
V in eq. ((S} (Bound approx.) and usdb, while by continuum 
states for V in eq. Q (CCSM). 
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sistently with experiments [20, |2l|] . Peak energies of vari- 
ous continuum spectra are depicted in Fig. [2l The CCSM 
results for ^'*0 agree very well to the experimental values 
reported recently by Hoffman et al. [7|. Note that the 
present calculation is done without adjustment to these 
experimental neutron spectra and also that it has been 
done before the experiment [2^ . 

We next evaluate the magnitude of the continuum- 
coupling effect. For this purpose, we first show the re- 
sults of a normal SM calculation performed by still using 
V in cq. (m but replacing the discrctizcd continuum ba- 
sis states, with the 0^3/2 HO wave function. In 
other words, this calculation is the same as the CCSM 
calculation except that there is no continuum coupling, 
and will be referred to as bound approximation. Vertical 
bars in Figure [1] stand for such results. 

Figure [T] (a,b) indicates that J=l and 2 continuum 
peaks appear, respectively, 1.4 and 1.2 MeV below the 
corresponding "bound approx" levels. This is purely due 
to the continuum coupling, as the interaction is the same 
between the CCSM and the bound approximation. One 
thus finds that the continuum coupling lowers energies by 
more than 1 MeV. We emphasize that a good agreement 
to experiments is owing to this size of the continuum 
coupling effect. We note also that the J'^=l+ and 2+ 
states remain unbound after including such significant 
effects of continuum coupling. 

The spacing between J'^=l+ and 2+ peaks (levels) be- 
comes smaller in the CCSM than in the bound approx- 
imation. This difference is consistent with experimental 
data as seen in Fig. [5] [2^. The two peaks with J^= 
1+ and 2+ are split due to the residual interaction V 
in eq. ([2]), which works attractively for J'^=2'^ but re- 
pulsively for J^=l+. The weaker splitting in the CCSM 
result is due to a weaker effect of V for the wave functions 
of CCSM as compared to HO wave functions. Namely, in 



the former, even inner part is more spread than the HO 
wave function, making the effect of V smaller. Note that 
the repulsive effect for J'^=l+ disappears if we use a zero- 
range interaction because of the exclusion of the state of 
total spin, 5=1. The penetrability should be different 
between J'^=2+ and J''=l+ because of their splitting, 
giving smaller width to J'^=2+ than to J'^=l^ . Thus, 
the widths also refiect properties of the residual interac- 
tion. This effect is not spherical and hence may not be 
taken into account by renormalizing the potential barrier 
height. 

Figure [1] (a,b) includes the results obtained by setting 
V = 0, i.e., no residual interaction. This is referred to 
as "no int." in the figure. In this case, the neutron 
spectrum becomes identical to the neutron emission from 
the doorway state formed by the 0c?3/2 HO wave function 
on top of the ^^O core. This is nothing but the pure 
single-particle picture and the J'" dependence disappears, 
as is clear in Fig. [TJ 

Figure [5] exhibits also the levels of the fully mixed sd- 
shell calculation by SDPF-M. First of all, these levels 
are pretty high in energy. We can confirm that these 
levels are nearly identical to the corresponding levels of 
the bound approximation with V of Eq. ^ in the filling 
scheme, reinforcing the validity of the usage of V in the 
present case. 

Figure [5] includes the levels of the normal sd-shell cal- 
culation by usd-b interaction |24| . where continuum ef- 



fects seem to be included into the effective interaction to 
a certain extent. 



V. WAVE FUNCTIONS OF CONTINUUM 
STATES IN 2*0 

We next discuss the density distribution of two valence 
neutrons with one in continuum ^3/2 orbits. The density 
p{r) is defined in general as the expectation value of the 
operator 

p(r) = iP''{r)iP{r) ^ ^ u*„^{r)un{r)al^an, (11) 

m,n 

where Umi'i') is the spatial wave function of the single- 
particle state with quantum number m, and (a) de- 
notes usual creation (annihilation) operator. We consider 
hereafter the density with the angular integration : 



p(r) = / dQpir). 



(12) 



This makes sense because one neutron occupies always 
one of the discretized c?3/2 orbits, while the other valence 
neutron is in the s orbit. This quantity is expressed as 

Pk{r) = {4\pir)\4) 



|0i,(r)p+ 5]c'(*)c|,.)0.L(O0.rf('^) 
Pis{r) + Pkd{r). 



(13) 
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FIG. 3: (color online) Densities (xr^) of ^3/2 for peaks of (a) 
1+ and (b) 2+ of ^*0, and (c) 3/2+ (c) of ^^O in the CCSM 
calculation, compared to those of harmonic oscillator (thin 
line). Solid lines indicate densities (arbitrary unit) of the 
corresponding single-particle "resonance" with phase shift of 
7r/2 with the same energy. 



Clearly, pkd represents the contribution from discretized 
continuum orbits to the fcth discrete eigenstate. 

Figure [3] (a) and (b) show, respectively, r'^pudi'r), for 
the appropriate cigenstates denoted by k which arc close 
to 1+ and 2+ peaks in Figs. [1] (a) and (b). Both show en- 
larged densities in and near the oxygen nucleus, while the 
enlargement differs between (a) and (b). The first node 
of the radial wave function is somewhat further away for 
2+ than for 1+. These density distributions are com- 
pared to the one calculated by harmonic oscillator 0^3/2- 
The large differences between the CCSM and HO results 
are seen, as expected. In addition, one sees notable dif- 
ferences between 1+ and 2"*", which are consequences of 
the different configuration mixing in continuum due to 
the nuclear force. 

We next discuss how the CCSM wave functions can 
be compared to scattering states, particularly, resonance 
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FIG. 4: (color online) Densities (xr^) of ^3/2 of ^^O for 
2+ peak state in the CCSM calculation, compared to those 
of single-particle resonance (solid line) wave functions which 
have the phase shift of 7r/2 with the same energy. 



ones. For this purpose, we consider a linear combina- 
tion of the incoming and outgoing states of with 
the phase shift 7r/2. As only asymptotic behavior is rele- 
vant here, wave functions of a free particle are considered. 
The energy is fixed to be the same as that of the emit- 
ted neutron from the peaks of Figs. [T] (a) and (b). This 
phase shift corresponds apparently to the single-particle 
resonance. We simply assume it without dealing with 
the potential. The resultant densities are shown also in 
Figs. [3] (a) and (b) as S.P. "resonance". Note that the 
S.P. "resonance" density has no normalization factor, and 
only its shape matters here. 

The asymptotic behavior of the density of the S.P. "res- 
onance" is very close to that of CCSM densities at larger 
r as seen in Fig. [3] (a). For the 2+ state in ^^O, the 
S.P. "resonance" density shows some deviation from the 
CCSM density (sec Fig.[3](b)). It becomes basically iden- 
tical to that of CCSM density at even larger distance as 
shown in Fig. |4j while this slower convergence for the 2+ 
state is of a certain interest. Thus, the CCSM calcula- 
tion always gives the correct wave length for the emitted 
neutron asymptotically, and furthermore the wave func- 
tion at the peak has the same phase shift to what we can 
expect from single-particle resonance at its energy. 



VI. RESONANCE IN NEUTRON SCATTERING 

In this section, we examine the similarity between the 
spectrum of the present approach and that of the neutron 
resonance on ^^O nucleus. We can calculate the phase 
shift of the neutron scattering for the potential of eq. ^ 
without the wall. Figure [5] (a) shows the (tangent of) 
phase shift thus calculated, which exhibits a clear pattern 
of the resonance. Figure E] (a) shows two spectra. One 
is of them (dashed line) is the spectrum of a neutron 
emitted from the HO doorway state built on the 
^^O core. This is a straightforward application of the 
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FIG. 5: (color online) (a) Phase shift for 3/2+ in ^^O. (b) 
Effective phase shift for 1+ and 2+ states in ^*0. 



calculations made for the ^''O case. The other spectrum 
(solid line) is the cross section calculated from the phase 
shift 6 J shown in Fig. [5] (a); 
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FIG. 6: (color online) (a) Spectra of neutron in continuum 
^3/2 in ^^O. Green dashed line denotes the CCSM spectrum 
of emitted neutron through the Harmonic Oscillator door- 
way state. Red solid line implies the cross section of the 
single-particle scattering by the same one-body potential as 
in the CCSM calculation, (b) Similar comparison between the 
doorway-state emission and the phase-shift approaches for the 
1+ and 2+ states in ^'^O. 



Air 

a J = -p-(2/ + l) siT?5j, 



(14) 



where k is the usual notation of the momentum and J 
(l) denotes the total (orbital) angular momentum of the 
channel. In this paper, 1=2, as only d wave is consid- 
ered for continuum. The similarity is remarkable between 
the two results, while the vertical scale is adjusted. In 
the resonance state, the incoming particle with the res- 
onance energy is considered to stay inside the imcleus 
for the resonance life time. If this meta-stable state in- 
side the nucleus has a wave function somewhat similar 
to the HO wave function, the outgoing particle follow a 
spectrum resembling that of the HO doorway state emis- 
sion. The similarity between the meta-stable wave func- 
tion and the HO wave function is expected naturally be- 
cause they reflect attractive mean effects of the nucleus, 
but these wavefunction cannot be identical either. The 
meta-stable one should be more diffuse. 

It is then of interest to see the cases of ^^O. Since 
the 1+ and 2+ states in ^^O are not of single-particle 
nature, we introduce effective phase shift first. The wave 
functions of two valence neutrons can be rewritten, by 



using eqs. (|7l8p . as 



where 



d3/2;J,k) =J2^i 1*^3/2) 



(15) 



The effective phase shifts are obtained from asymptotic 
behavior of the wave function of the state 1^3/2; j.fe)- Fig- 
ure [S](b) shows the (tangent of) effective phase shift thus 
obtained. The effective phase shift behaves quite simi- 
larly to the phase shift of single-particle state, while it 
makes sense only in a particular channel, which is speci- 
fied by J in the present study. The actual wave function 
of a particle in continuum corresponds to a certain su- 
perposition of wave functions of different channels, and 
therefore there is no unique phase shift describing the 
scattering. We can still use the effective phase shift for a 
particular channel very conveniently. 

The "resonance" energy can be shifted due to the cor- 
relation by the interaction in eq. Figure EKb) shows 
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four spectra. Two of them (dashed hne) are those shown 
in FigdJ The other two (soUd hne) are the cross sec- 
tions calculated by eq. (fT^ with the effective phase shifts. 
The spectra obtained by the two methods appear to be 
similar to each other, while certain differences are visi- 
ble also. We stress that these two methods are for dif- 
ferent processes and certain differences are rather nat- 
ural. Table |T] summarizes peak energies and widths of 
the doorway-state emission spectra (shown as "CCSM" ) 
and the phase-shift spectra (shown as "phase shift"), in 
comparison to experiments. 

The width of the doorway-state emission spectrum 
shows larger widths than the corresponding phase-shift 
spectrum (see also Fig. [5]). The difference here is con- 
sidered to stem from the difference between the doorway 
state (HO wave function) of the knockout reaction and 
the inner wave function of the phase-shift approaches. In 
the former, radial wave functions are more confined inside 
the nucleus, making the mixing with various continuum 
states stronger and the width larger. Table |T] indicates 
that the difference of the width is largest for the 2+ state 
in ^^O. This fact is consistent with the above argument as 
the 2+ state has the lowest peak among the peaks shown 
m Table HI and its inner wave function should be most 
different from HO wave function. 



VII. LOW- LYING STATES IN ^^O AND ^^O 



We now discuss the continuum properties of ^^O and 
^^O. As for ^^O, there are three neutrons in the present 
model space. We diagonalize the Hamiltonian in eq. 



using the basis, |(lsi/2)^ (8) 1^3/2; 3/2+) (z = 1, ■ 



Namely, one neutron is placed in the continuum ^3/2 
space, for simplicity. We discuss the reaction of single 
proton knockout from ^^F, namely, ^Be(^^F,^^ 0)X [3]. 
Here we assume, by a similar argument to the ^^O case, 
that the neutron emission probability is proportional to 

p)^^ calculated in this configuration space. The resultant 
peak energy is included in Fig. [21 exhibiting a good agree- 
ment to experiment 0]. The spectrum of the emitted 
neutron is shown in Fig. [TJc) . The density distribution 
of this continuum state is shown in Fig. |3l 

The ground and the first 2+ states of ^^O are calcu- 
lated within the basis states: |(lsi/2)^ ® ^^3/2' (* ~ 
Ij ■ ■ ■ jimax)- These states turn out to be unbound in 
agreement with experiment [lol . [ll| . The resultant peak 
energies are included in Fig.|51 The 0+-2+ spacing in the 
CCSM result is smaller than in the normal SM, which 
may indicate a weakening of the pairing gap. Note that 
the density dependence of the pairing interaction is not 
considered here. The results obtained by the usdb inter- 
action [l^l are shown too. Note that a part of continuum 
effects are included in the usdb TBME's by fit. 

Table |T] includes the calculated results and the experi- 
mental data of 250 and ^^O 0, [H, [2|. 



TABLE I: Peak energies and widths of spectra of emitted neu- 
trons. Peak energies are measured from the particle-emission 
thresholds [MeV]. "CCSM" refers to spectra of emitted neu- 
tron from HO doorway states, while "phase shift" to those 
of neutron scattering calculated in the phase-shift approach. 
Results of the normal shell model calculation in the full sd- 
shell and the bound-state approximation are shown. GSM 
results with the same Hamiltonian are exhibited. 
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phase shift E 
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FIG. 7: (color online) Convergence of the results with respect 
to the position of the infinite wall, L. 



VIII. CONVERGENCE OF CCSM WALL 
POSITION 



We shall comment on the convergence of the results 
with respect to the position of the infinite wall, L. In 
Fig. [3 the peak energies discussed already for oxygen 
isotopes are plotted as a function of L. It is seen that 
the results do not change so much if the value of L is 
taken to be sufficiently large. It should be noted that 
the results with even usual values of L ~ 50 fm are not 
stable. Note that the value of L is fixed to 3000 fm in 
this paper. 
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IX. GSM CALCULATION BY THE SAME 
HAMILTONIAN 

The present continuum properties can be studied by 
other methods, e.g., Gamow Shell Model (GSM) 
and Coupled Channel calculation [2^. The shell model 
in continuum (CSM) [29| is another approach where a 
non-hermitian energy-dependent Hamiltonian is used, by 
evaluating continuum coupling through one-body poten- 
tial. 

We here investigate the relation of the GSM [2^,[13] to 
the CCSM using the present Hamiltonian, Eq. ([3]). 

The GSM is one of the extensions of the normal shell 
model to continuum, by introducing the so-called Gamow 
basis. The Gamov basis is generated by analytically con- 
tinued Schrodinger equation 



S.P. resonance 

free particle with phase shift nJ2 




r[fm] 

FIG. 9: (color online) Single-particle resonance state for ^^O 
and a free-particle state with 5i = 7r/2. 



c+ 



dqq^Vi{k,q)%l)ni{q) = E^i-ipnlik) (16) 



where fi is the reduced mass between a neutron and ^^O 
core in this study. VJ(fc,fc') is the Fourier-Bessel trans- 
form of spherical Woods-Saxon potential which is the 
same as in Eq. ([T]) except for the wall. 



Re(k)fm-1 




Re (k)fm-l 



contour L 



2 



FIG. 
basis. 



8: Two contours used for the construction of Gamow 



TBME's of the effective interaction in Eq. ((2| are cal- 
culated following the prescription in Diagonalizing 
the Hamiltonian matrix which is complex and symmetric, 
we obtain the cigenstate whose imaginary part is inter- 
preted as the width of the resonance state. Note that the 
GSM provides us with resonance properties. This calcu- 
lation produces practically the same peak energies but 
smaller (by 5%^30%) widths (see Table |l| as compared 
to the CCSM calculations with the same Hamiltonian. 
The CCSM phase-shift calculations produce results in 
between or closer to CCSM doorway-state-emission cal- 
culation. 

Since we discretizc the integral equation ()16p . we have 
to check the convergence of the GSM results as a function 
of the number of single-particle states belonging to the 
Gamow basis. We confirmed that converged results can 
be obtained with 30 basis states. 

In the GSM calculation, there is a degree of freedom to 
choose the deformed contour. We use two different types 
of the contours as shown in Fig. |8l specified by the points 
A{0.3 - 0.2i fm"^), 5(0.5 fm"^) and C(0.3 - 0.2ifm"^). 
It has been confirmed that the result does not change 
much between these two contours. 



X. METHODOLOGICAL RELATION BETWEEN 
CCSM AND GSM 

We shall look into the methodological relation of the 
CCSM calculation to the GSM. 

The resonance energy can be determined in general 
by two different ways. One is to determine it as a point 
where the phase shift goes through 7r/2, as has been done 
with CCSM wave functions. The other one is to treat a 
resonance state as a pole of the S-matrix. The latter 
definition is equivalent to the situation that the radial 
wave function of the resonant state has the pure outgo- 
ing boundary condition [33| . These two are of course 
related to each other, but their mathematical relation, 
which may have been noticed, appears to be not so sim- 
ple, as discussed below. The asymptotic behavior of the 
radial wave function is given by the combination of the 
Hankel functions 



Ri{r) Ah\^\kr) + Bh\^\kr), 



(17) 



where k=y/2fj,E/h. The coefficients can be expressed as 
^^g2j5,(fc)y2 ^ Si{k)/2 and B=l/2. One then obtains at 
large r limit. 



Ri{r) e'^' [cosSiji{kr) + sm5ini{kr)] 



(18) 



At the resonance state defined by (5; = 7r/2, i?;(r) is rep- 
resented as a standing wave, but not an outgoing wave as 
in the case of complex approaches. The asymptotic wave 
functions discussed in Figs. |3] and |4] correspond to this 
situation. Figure |5] shows the relation between the wave 
function of the single-particle resonance of ^'^O shown in 
FigE] (a) and that of the free particle of the phase shift 
Si = Tr/2. 

If the phase shift Si sharply increases by about tt in the 
vicinity of A: = ko{k,ko S R), going through 7r/2 at fco, 
the energy dependence may be approximated by 



tan S] 



ki 



ko — k ' 



fci > 0. 



(19) 
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It then follows that near k = kg, 



sin' 6i = ^\Si 



kl 



rV4 



{k - koy 



kl 



(20) 



where Eq and T correspond to the peak energy and the 
width of the peak, respectively. Combining with eq. ([T^ , 
this leads us to the Breit-Wigner formula, with 



Ec, = 



2^' 



r = 



2koki 



(21) 



Another usua l ap proach of the resonance phenomenon 
is the S-matrix [3ll]. The S-matrix takes, for real fc's, the 
form 



Si{k) 



k — ko — iki 

k — ko + iki 



(22) 



When this S-matrix is analytically continued to the com- 
plex momentum-plane. Si has a pole at fc = fco — iki in 
the lower half of the complex k-plane. We then obtain Eq 
and r, as discussed in |34| . from the pole of the S-matrix, 



En 



^e^(fco - ifci)^ 



(/cq kl) 



r = — 2/m — (fco — iki) 
2fi 



1 

2 2fcofcl 



(23) 
(24) 



Thus, the pole of S-matrix corresponds to the resonance 
state, only if this pole in the complex k-plane is close 
enough to the positive real axis and well isolated from 
other poles. If one applies the width for the 1+ state 
of ^^O shown in Table U the estimated shift of the peak 
energy turn out to be about 0.004 MeV, which is small 
enough. 

The peak energies and widths obtained by the CCSM 
doorway-state-emission and phase-shift calculations can 
be compared with GSM results in Table U for the same 
Hamiltonian. At first glance, the peak energies are in 
a good agreement among the three calculations. The 
width becomes smaller in going from the CCSM doorway- 
state-emission to the CCSM phase-shift, and to the GSM 
calculations in general. The difference may be partly due 
to the coupling/mixing between basis states evaluated 
differently among the three, and also due to the different 
boundary conditions as discussed. 

We comment on the relation to Ref. [s^, where the 
microscopically derived two-body interaction is used for 
GSM calculation. The oxygen isotopes were studied. 
As well known, the direct use of microscopic interaction 
into the shell model calculation does not work for bound 
states. The same feature is seen in continuum calcula- 
tions in Ref. [1^. This problem has been avoided in the 
present work, focusing on the continuum coupling effect. 
We note that the present work has been started in 2002 



[371 ] , and has been refined in 2008 [2^ . The present work 
does not owe to results of [1^1, as it has been completed 
before 

XI. SUMMARY 



In summary, we introduced the Continuum-Coupled 
Shell Model (CCSM). The CCSM Hamiltonian is de- 
signed to include correlations due to a two-body nuclear 
force in the continuum, in a manner consistent with the 
realistic SM interaction for bound states. Spectra of 
emitted neutrons from knockout reactions are calculated, 
in good agreements to recent experiments without adjust- 
ment to final results. The continuum effects are signifi- 
cant; the difference from bound-state calculation is more 
than 1 McV for the cases studied. We point out that the 
1+-2"*" spacing of ^'^O is smaller in the CCSM than in the 
bound-state calculation, and that the CCSM result in- 
deed agrees well with the recent experiment (3] . This sug- 
gests that the appropriate treatment of the nuclear force 
is still important in the continuum. We also calculated 
the resonance states by a phase shift approach and the 
Gamow Shell Model (GSM) using the same Hamiltonian 
(Eq. Q). The peak energies agree with those of CCSM 
doorway-state emission calculation rather will, although 
the doorway-state emission and the neutron-scattering 
resonance are very different processes in principle. The 
width of the resonance states were obtained to be smaller 
than the case of CCSM doorway-state-emission calcula- 
tion. This difference originates partly in the differences 
between the doorway and the resonance states. On the 
other hand, the differences among the CCSM (doorway), 
CCSM (phase shift) and GSM appear rather minor, par- 
ticularly as far as the peak energies are concerned. 

The present work suggests that the high single-particle 
energies for the bound 0^3/2, which are used for shell- 
model calculations for bound states, are indeed consistent 
with observed continuum properties. If they were not 
high enough, the continuum spectra would become too 
low compared to experiments. There should be a strong 
shell evolution mechanism so that this orbit comes down 
in F isotopes [1, H, |3^ . As the CCSM idea naturally ex- 
tends the successful SM interaction to continuum states 
(or open quantum systems) , its future developments into 
full configurations and more applications to exotic nu- 
clei are intriguing issues. For instance, the continuum 
"ground state" of 26,28q attracts much interest as a pos- 
sible "mother of pearl" of neutron nugget /matter. The 
continuum effect was not taken into account for bound 
states in this paper for simplicity. There may be some 
cases where the continuum couplings can be crucial even 
for bound states, e.g., ^^Li. 

We thank Prof. K. Kato for valuable comments. This 
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